Hot QCD equations of state and relativistic heavy ion collisions 
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We study two recently proposed equations of state (EOS) which are obtained from high temper- 
ature QCD, and show how they can be adapted to use them for making predictions for relativistic 
heavy ion collisions. The method involves extracting equilibrium distribution functions for quarks 
and gluons from the EOS, which in turn will allow a determination of the transport and other 
bulk properties of the quark gluon plasma. Simultaneously, the method also yields a quasi particle 
description of interacting quarks and gluons. The first EOS is perturbative in the QCD coupling 
constant and has contributions of 0[g''). The second EOS is an improvement over the first, with 
contributions upto 0{g^ln{^)); it incorporates the nonperturbative hard thermal contributions. The 
interaction effects are shown to be captured entirely by the effective chemical potentials for the glu- 
ons and the quarks, in both the cases. The chemical potential is seen to be highly sensitive to the 
EOS. As an application, we determine the screening lengths which are, indeed the most important 
diagnostics for QGP. The screening lengths are seen to behave drastically differently depending on 
the EOS considered., and yield, therefore, a way to distinguish the two equations of state in heavy 
ion collisions. 
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I. INTRODUCTION 

Recent experimental results [H H, H, 0] indicate that 
the quark gluon plasma has already been produced at 
RHIC, and that its behavior is not close to that of an 
ideal gas. Indeed, measurements of flow parameters [l| , 
and observations of jet quenching[5j] have stimulated the 
theoretical interpretation that the QGP behaves like a 
nearly perfect fluid characterized by a small value of 
the viscosity to entropy density ratio, lying in the range 
.1 ~ .3. [3, H, Q; this range may be contrasted with 
the corresponding value for liquid Helium (above super- 
fluid transition temperature) which is close to ten 'lO|. 
These observations signal the fact that the deconfined 
phase is strongly interacting, and are consistent with the 
lattice simulations 11| , which predict a strongly inter- 
acting behavior even at temperatures which are a few 
Tc- In an attempt to appreciate this surprising result, 
interesting analogies have been drawn with ADS/CFT 
correspondence ll(| and also with some strongly coupled 
classical systems 12]. In any case, the emergence of the 
strongly interacting behaviour puts into doubt the credi- 
bility of a large body of analyses which are based on ideal 
or nearly ideal behaviour of QGP. 

In this context, there is an interesting att emp t by 
Arnold and Zhai 13] and Zhai and Kastening IJ] who 



have determined the equation of state (EOS) of interact- 
ing quarks and gluons upto 0{g^) in the coupling con- 
stant . This strictly perturbative EOS, which we hence- 
forth denote by EOSl, has been improved upon by Ka- 
jantie et al.[l3[3 '^^'^ have incorporated the contribu- 
tions from the nonperturbative scales viz. qT and g^T 
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and determined the EOS upto 0{g^ln{^))^. The latter 
will be denoted by E0S2. Subsequent studies [ii,[23,[2ll 
have emphasized the relevance of the above EOS to study 
quark gluon plasma. One would naturally wish to com- 
pare these EOS with (the fully non-perturbative) lattice 
results. E0S2 has been found to be in qualitative 
agreement with the lattice results. It is not without in- 
terest to explore further whether this qualitative agree- 
ment can be further quantified, and whether HTL im- 
proved EOS can describe the QGP produced in heavy 
ion collisions. It is worthwhile noting the earlier attempts 
[23 , [23l . [2^ that have been made to determine thermo- 
dynamic quantities such as entropy and the specific heat 
Cv in improved perturbative approaches to QGP. 

On the other hand, it is by now well established that 
the semiclassical approach is a convenient way to study 
the bulk properties of QGP [H [H, HHH, El , since they 
automatically incorporate the HTL effectspsl. I29I ] . There 
is a wealth of results which have been obtained within 
this framework j25l [26l . [27j . where the nonperturbative 
features manifest as effective mean color fields. These 
color fields have the dual role of producing the soft and 
semisoft partons, apart from modulating their interac- 
tions. The emergence of such effective field degrees of 
freedom, together with a classical transport has been in- 
dicated earlier by Blaizot and lancu [30j . 

In this context, it is pertinent to ask if one could use 
heavy ion collisions to distinguish the various EOS and 
pick up the right one, by employing the semiclassical 
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framework involving an appropriate kinetic equation. The 
purpose of this paper is to explore such possibilities. As 
a first step in this direction, we shall show how the dis- 
tribution functions underlying the proposed EOS can be 
extracted with a minimal ansatz, viz, effective chemical 
potentials for quarks and gluons. Once the distribution 
function is obtained, it can be used to study the bulk 
properties of the system such as chromo responses in- 
cluding the ubiquitous Debye Mass. Postponing all the 
other applications to a future work, we shall concentrate 
on determining the Debye mass through this procedure. 
As mentioned, we focus on EOSl and E0S2. Both of 
them have been proposed for the case when the baryon 
number density vanishes. The corresponding chemical 
potentials are hence set to zero. There exist generaliza- 
tions of the above EOS, proposed by Vuorinen and 
more recently by Ipp et al [18], which allow for a finite 
baryon number. The two sets arc applicable to distinct 
physical situations; the former (EOSl and E0S2) are 
relevant to the QGP in the central midrapidity region of 
URHIC while the works of Ref. [13, Ell are applicable to 
peripheral collisions and /or when the so called nuclear 
transparency is only partial. An application of the above 



EOS to URHIC will be taken up separately. 

This paper is organized as follows. In the next section, 
we extract the distribution functions for the gluons and 
the quarks from EOSl and E0S2. Wc consider the pure 
gluonic case separately from the full QCD, by first set- 
ting Nf = 0. The (interacting) quark sector is then dealt 
with. In section III, the Debye mass is determined by 
employing the semiclassical method developed by Kelly 
et.al [2^, l29|. In section III(B), we compare our results 
on screening length for EOSl and E0S2 with the recent 
lattice results. We summarize the results and conclude 
in section IV. The appendix contains some details of cal- 
culations which are not explicitly given in the main text; 
it also lists some useful integrals. 



II. EXTRACTION OF THE DISTRIBUTION 
FUNCTIONS 

Recently Arnold et. al have derived an equation 
of state (EOSl) for high temperature QCD up to 0{g^). 
EOSl reads. 
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EOSl has been subsequently improved by Kajantie et.al. 
[isl \ld\ who proposed another equation of state (E0S2) 
by improving the accuracy to the next order in the cou- 
pling constant, and also included the HTL effects; recall 



that the latter are essentially nonperturbative, and con- 
tain contributions from scales T, gT, and g^T. E0S2, 
which is thus determined upto 0{g^ln-), has the form 
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In the above expressions, Nf is the number of fermions, the renormalization scale parameter in the MS scheme. 
Us = (7^/(47r) is the strong coupling constant and /xms is 
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Note that as runs with /3 and flms- As remarked, the 
utihty of this EOS in the context of QGP thermodynam- 
ics has been discussed ear her by Rebhanplj. 

We now set to determine equihbrium distribution func- 
tions < Hgj > for the gluons and the quarks such that 
they would yield the EOS given above. The ansatz for 
the determination involves retaining the ideal distribu- 
tion forms, with the chemical potentials fig and fif being 
free parameters. Note that for the massless quarks ( u 
and d) which we consider to constitute the bulk of the 
plasma, fi = if they were not interacting. This approach 
is of course not novel, since it underlies many of the ideas 
that attempt to describe the interaction effects in terms 
of the quasiparticle degrees of freedom. In the present 
context, we refer the reader to Ref. [sil, [s^l, where an 
attempt is made to describe the lattice results in terms 
of effective mass for the partons. 

We pause to note that the chemical potentials which we 
introduce are not the same as those which yield a nonzero 
baryon number density, as e.g., in [13, Here, the 

chemical poetntial merely serves to map the interacting 
quarks and gluons at zero baryon number chemical po- 
tential to noninteracting quasiparicles, viz., the dressed 
quarks and gluons. Their interpretation is, therefore, 
more akin to the effective mass, albeit as functions of 
the renormalization scale and temperature as we show 
below. Thus, the baryon number density of the plasma 
continues to vanish. 

As the first step in our approach, we express the EOS 
in the form 

P^P^ + P^ + APg + APf (3) 



The first two terms in the RHS of Eq.Q are identified 
with the distributions of an ideal gas of quarks and glu- 
ons. The effects of the interaction in pure QCD are rep- 
resented by APg and the residual interaction effects, by 
APf. For the EOS which we are interested in, the identi- 
fication of the above terms is straight forward. APg can 
be identified by first setting Nf = and then subtracting 
the ideal part. The residual term is naturally identified 
as APf after subtracting the ideal part for quarks. In 
general the form of EOS (see Eq.([l]) and Eq.©) 
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For EOSl the coefficients A, B, C, D are denoted with 
a prime and are given by 
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whereas for E0S2 the coefficients can be written in terms of the above primed coefficients for EOSl, as 

AiasiNf)) = A'(a,(iV,)) + 1134.8(^)3 log(-) 

TT as 

B{as{Nf)) = B'{asiNj))-^i^)Hog{-) 

Z TT as 

Cia.m.N,) = C-,a,,iV,).«,) + (65.89iV, + 7.653A'J)(^li)Mog(-L) 



4 



D{a,{Nf),Nf) = D'ia,{Nf),Nf) 



1485,, 1 2 , -..a^.o 1 , 

^[(1 + -iV/)(l - -Nf) - l](-^)3log(-). 
ZD oo n as 



(7) 



We seek to parametrize the contributions from all the 
non-ideal coefficients in terms the chemical potentials /ig 
and fif for gluons and quarks respectively. Since the EOS 
have been proposed at high T, with their validity being 
at temperatures greater than 2Tc [33| . we treat the di- 
mensionless quantity flgj = PfJ-gj perturbatively. This 
approximation needs to be implemented self consistently, 
and accordingly, we expand the the grand canonical par- 
tition functions for gluons and quarks as a Taylor series 
in fig J- We obtain the following expressions: 
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where, Zg and Zf are given by 
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We determine Zg and Zq (defined in Eq.®) up to 
0{jlgjY. The truncation is seen to yield an accuracy 
of ~ ten percent when we consider EOSl. However, for 
the more physical E0S2, the agreement is within one per- 
cent. The gluon chemical potential /ig gets determined 
by 
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Similarly, the equation determining reads: 
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The coefficients Ag"-* and A^^^ are given by 

4")=9^^1og(Z,)|(,^=o) 
4")=atyog(Z/)|(,,^.o) (12) 

The explicit forms of A" up to n = 3 are listed in the 
Appendix. 



A. Explicit evaluation of the chemical potential 

Before we discuss the solution of Eq. ([TU]) and Eq. pTjl , 
it is instructive to evaluate ^gj with just the linear and 
quadratic terms for the sake of comparison. In the linear 
order, the solutions read 
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while, in the next to the leading order, the solution has 
the form 
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where, Cg = 2APg(3V / A'g^'^ and C/ = 2APfPV/Af\ 
Finally, the exact solutions can be obtained by using well 
known algebraic techniques. Since the explicit algebraic 
solutions do not have an illuminating form, we show the 
solutions graphically instead in the next subsection. 

The distribution functions for the gluons and the 
quarks get determined, in terms of the chemical poten- 
tials, through 
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The extraction of the distribution functions is, never- 
theless, incomplete. For, the EOS -and hence the chemi- 
cal potentials - depend on the renormalization scale. On 
the other hand, the physical observables should be scale 
independent. We circumvent the problem by trading off 
the dependence on /Ims to a dependence on the critical 
temperature Tc. To that end, we exploit the temperature 
dependence of the coupling constant (T) [3I, [11] and 
of renormalization scale: 
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where, bo = 33-2Nf/12n and Xms = l-14Tc. With 
this step, the distribution functions get determined com- 
pletely, and are obtained as functions of T/Tc. 

We note that the results presented below, being valid 
for T > 2Tc, need to be supplemented by a similar analy- 
sis for EOS which are valid for T ^ Tc- Such an analysis 
does indeed exist, along the lines of this paper (3^ . who 
have considered the lattice EOS. They do not determine 
the Debye mass, but focus on the impact of the EOS on 
the flow parameters in heavy ion collisions. 

Although, EOSl and E0S2 have been computed 
within the framework of weak coupling technique but 
they give convergent results for the temperature ranges 
which are more than 5Tc. We shall see in the next Sub- 
section that these equations of state are far away from 
their ideal behaviour up to the extent that they can be 
utilized to make definite predictions for QGP. 

B. Hot QCD EOS vs Ideal EOS 

As a warm up, we compare EOSl and E0S2 with the 
ideal EOS by plotting the ratio R = pi^pi , as functions 

of temperature, in Figs. 1,2. 
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FIG. 1: (color online) Behaviour of R with temperature for 
EOSl 

The most striking feature that we see is the large sen- 
sitivity to the inclusion of the g^ln{^) contributions. It 
is most pronounced in the behaviour of pure QCD where 
major qualitative and quantitative differences appear: (i) 
For EOSl, R increases with T, in contrast to E0S2 which 
where it decreases from above, approaching the same 
asymptotic value for large T; (ii) interestingly, the non- 
perturbative (and higher order) corrections makes the 
system less non-ideal. Indeed, EOSl yields values of R 
which are 10 — 45 percent away from the ideal value 1, in 
contrast to E0S2 (see Fig2), for which R is only 2 — 8 per- 
cent away from the ideal value. Incidentally, the above 
observation implies that the expansion in Eq.(I5]) works 
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FIG. 2: (color online) Behaviour of R with temperature for 
E0S2 



better for E0S2 than the former. This will be reflected 
later in the behavior of effective chemical potentials with 
temperature. We shall further see that, smaller the value 
of l/igji, better will be the approximation. 



C. The chemical potentials 

The variation of the effective chemical potentails with 
renormalization scale at a fixed temperature has already 
been studied in reference Il4| . We prefer to recast it into 
a dependence of p,gj on ^ (see the previous subsection) 
since it is more relevant to the study of QGP in heavy 
ion collisions. This is shown in Figs. 3 — 8, where the 
contributions coming from linear, quadratic and cubic 
approximations in the Taylor series(Eq.([5])) are individ- 
ually displayed, for both EOSl and E0S2. These figures 
exhibit, in essence, all the interaction effects. 

A number of features emerge from examining Figs. 
3 — 8. Consider EOSl first. Here, the linear approxima- 
tion does reasonably well for pure QCD, but fails badly 
in the quark sector. The chemical potential is negative in 
both the sectors, and approaches the ideal value asymp- 
totically from below. In contrast, E0S2 leads to a dif- 
ferent behaviour: fig starts with a small positive value 
at T ~ 2Tc, and stays essentially so until T ISTc, and 
switches sign to acquire a small negative value. Since the 
magnitude remains less than 0.1 throughout, the devia- 
tion from the ideal behaviour is minimal, fl / remains neg- 
ative (with the maximum magnitude ~ 0.25 at T = 2Tc), 
which is about a factor four smaller in comparison with 
the corresponding value from EOSl. The interaction ef- 
fects get manifestly stronger as we increase the number 
of flavors. 

It is significant that the ideal value is not reached even 
at T ~ lOTc, which indicates that the phase remains 
interacting. We also note that our method of extracting 
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FIG. 3: (color online) Effective chemical potentials at the 
linear order for EOSl 



FIG. 5: (color online) Effective chemical potentials at the 
quadratic order for EOSl 
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FIG. 4: (color online) Effective chemical potentials at the 
linear order for E0S2 
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FIG. 6: (color online) Effective chemical potentials at the 
quadratic order for E0S2 



the chemical potential works more efficiently for E0S2, 
as indicated by small corrections from higher order terms 
to the linear estimate. 



III. THE DEBYE MASS AND SCREENING 
LENGTH 

The extraction of the equilibrium distribution func- 
tions affords a determination of the Debye mass, via the 
semiclassical transport theory [1^ . The Debye mass con- 
trols the number of bound states in heavy qq systems, 
yields the extent of J/^ suppression in heavy ion colli- 
sions, provided that we have a reliable estimate of the 
temperature of the plasma. Even otherwise, the qual- 
itative significance of the Debye mass cannot be over 
estimated since the deconfined phase remains strongly 



interacting even at large T. 

The determination of mn is straightforward if we em- 
ploy the classical transport theory [29] . It is simply given 

by 

= g'^Cgj / ^ < ngj > d^p. (19) 

The above expression has to be used cautiously, though. 
The coupling constant g' in eq ll9l has a phenomeno- 
logical character, and should not be confused with the 
fundamental constant g appearing in the EOS. Keeping 
this in mind, we recall that if the plasma were to be 
comprised of ideal massless partons, the Debye mass 
would be given by 
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is obtained as 
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FIG. 7: (color online) Effective chemical potentials at the 
cubic order for EOSl 
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The hot QCD EOS modify the above expression; It 
is easy to see, from Eqs. dT^ that the new Debye 
masses, scaled with respect to their respective ideal 
values get determined, in terms of the standard PolyLog 
functions [43II by 
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It is, however, more convenient to plot the inverse debye 
mass i.e, the screening length as a function oiT/Tc- 



A. Relative screening lengths 

We first establish the notations. Let Xh denote the 
screening length generated by the hot EOS. Let A^^ be 
the screening length of an ideal qgp. It is convenient 
to consider also the contribution coming from the pure 
QCD sector, whose screening lengths we denote by A)^ 
and Af^ respectively. 
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Consequently, the expression for the total relative mass 



FIG. 9: (color online) The relative Debye Screening Length 
Tlh/g for EOSl as a function of temperature. Note that it is 
> 0.7. 



The behaviour of the screening lengths is shown in 
Figs. 9 -12. As in the case of the chemical potentials, the 
dependence on the order of perturbation is striking here 
as well. For EOSl where the contributions upto 0{g^) 
are included, the screening lengths in the full QCD as well 
as pure QCD remain nonzero. The dominant contribu- 
tion is from the gluonic sector, which dominates over the 
quark sector, as may be seen in Fig. 9 where we plot the 
ratio Ti-h/g — A;i/A^, which is in excess of 0.7 through- 
out. Note, however, that the relative dominance gets 
weaker as we increase the number of flavours. Fig. 10 
shows the variation of the ratio TZh/id = Xh/Xid as a 
function of temperature. Interestingly, the interaction is 
seen to weaken the screening, and so does an increase in 
the number of flavors. 

These results are in sharp contrast with the case of 
E0S2, which we recall has nonperturbative 0{g^ln{^)) 
contributions. These are shown in Figs. 11 and 12. It is 
clear from Fig 11 that the contribution from the pure glu- 
onic sector saturates the contribution to the screening all 
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FIG. 10; (color online) The relative screening length7?,(,/id for 
EOSl as a function of temperature. 



FIG. 12: (color online) The relative screening length TZu/id 
for EOS2 as a function of temperature. 
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FIG. 11: (color online) The relative Debye Screening Length 
TZh/g for E0S2 as a function of temperature. Note that it 
stays at 1 all the way upto ISTc. 



the way upto temperatures T ~ 13Tc, and drops sharply 
thereof. This feature is reinforced by Fig. 12 where the 
ratio Ti-h/id stays at zero between 2Tc and 12 — 13Tc. It 
is of a purely academic interest that the screening length 
should become non-zero beyond 12Tc. 

It appears that the perfect screening is indeed the 
strongest prediction of E0S2 and must be most easily 
tested in heavy ion collisions, where temperatures upto 
3Tc are expected at LHC. This is in sharp contrast with 
the assumptions of a near ideal behaviour, and also some 
theoretical analyses which in fact propose an enhanced 
production of J/^E* at LHC energies [40|. We, therefore, 
attempt to compare these predictions with the lattice re- 
sults below. 



B. Comparison with the lattice results 

In this subsection, we compare our results on screening 
length of EOSl and E0S2 with the lattice results. Lat- 
tice compuattions extract the screening lengths from the 
quark-antiquark free energies. To be concrete, we make 
the comparison with three distinct values of the coupling 
constant, g' = .3, .5, .8. 
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1 
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full QCD with Nf=2 
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FIG. 13: (color online) Behaviour of 2/Md with T/Tc for 
g' — .3 for EOSl. Note that 2/Mo is measured in fm 



Further, we consider three cases, (i) pure QCD, (ii) 
Np — 2 and (iii) Np = 3. To facilitate a proper com- 
parison, we take the respective transition temperatures 
to be Tc = 270MeV, 203MeV and 195MeV, as given 
by lattice computations. The comparison is shown only 
with EOSl since E0S2 predicts absolute screening in the 
range 2Tc < T < \2Tc that we are interested in. The re- 
sults are shown in Figs. 13-15. As observed, the screen- 
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FIG. 14: (color online) Behaviour of 2/A/15 with T/Tc for 
g' — .5 for EOSl. Note that 2/Md is measured in fm 



FIG. 15: (color online) Behaviour of 2/Md with T/Tc for 
g' — .8 for EOSl. Note that 2/Md is measured in fm 



ing weakens with increasing g'; moreover, the screening 
weakens with the increase in the number of flavours as 
well. For an explicit comparison, we consider the results 
reported by Kaczmarek and Zantow[4ll| who detemine 
the screening length by identifying it essentially with the 
first moment of the gq free energy. Their results are dis- 
played in Fig. 2 of [4l| to which we refer henceforth. 
Interestingly, the same qualitative features are exhibited 
EOSl and the lattice results, in both the aspects, viz, 
the dependence on the coupling constant as well as on 
the number of flavours. However, the agreement fails 
to get quantitative. The lattice results predict screening 
lengths which are smaller in value, except for Np = 3 
than the EOSl results. Indeed, the lattice screening 
length is ^ 0.7 fm in the vicinity of Tc, and drops to 
^ 0.4/m close to 2Tc. It is evident from Figs. 13-15 that 
the results of EOSl are 3 — 10 times higher in value. Any 
better agreement with a further increase in the value of 
g' is ruled out sinceg' < 1 necessarily. 



IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we have extracted the distribution func- 
tions for gluons and quarks from two equations of state, 
in terms of effective chemical potentials for the partons. 
The chemical potentials are shown to be highly sensitive 
to the inclusion of 0{g^ln{l/ g)) contributions, and exhib- 
ited most vividly by the screening length. Surprisingly, 
E0S2 which has interactions upto 0{g^ln{^)) shows less 
nonideal behaviour compared to EOSl (which has con- 
tributions upto 0{g^). Equally strikingly, the plasma 
corresponding to E0S2 is predominantly gluonic, in the 
sense that the Debye mass from the gluonic sector di- 
verges in the range 2Tc < T < \2Tc- This result is in 
contrast with the less precise EOSl, where the gluonic 
contribution is not that overwhelming. 



To place our analysis in perspective, we note that our 
analysis is based on but two equations of state, neither 
of which has full non-perturbtive contributions. Nev- 
ertheless, it may not be without merit since E0S2, for 
instance, makes rather strong predictions which may be 
tested; EOSl is, on the other hand, seen to be in qualita- 
tive agreement with the lattice results. Indeed, the work 
does provide a platform to study quantitatively the im- 
port of the EOS to heavy ion collisions in a quantita- 
tive manner. Experiments at LHC may be able to probe 
these EOS since a temperature in the range T ~ 2 — STc 
is expected to be achieved there. More importantly, the 
method developed here can be easily employed to study 
more precise EOS (as from lattice computations), or more 
general EOS (as the inclusion of baryonic chemical po- 
tentials). To be sure, an incisive analysis is possible 
only after studying other quantities such as the viscosity, 
its anomalous component [H, [s^, the viscosity to en- 
tropy ratio, and the specific heat. Finally,the insertion of 
the appropriate equilibrium distribution functions in the 
semiclassical transport equations allow for studying (i) 
the production and th e eq uilibration rates for the QGP 
in heavy ion collisions f25j, [2^ [27| , and (ii) the color re- 
sponse functions [42]. of which the Debye mass is but one 
limiting parameter. These will be taken up in subsequent 
publications. These investigations will be taken up sep- 
arately. 
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VI. APPENDIX 



We use the following standard integrals while extract- 
ing effective chemical potential: 



r^^TT^^^^T^^^ - 2(i+f:(^n^^±A^)) 

Jo (1 - exp(-p))3 (n + 1)3 J-^A fc! 

' ^ \ir^'n^3 -^^- ^+log(2) (23) 
(l + exp(p))3 12 

= 2 cxp(p) PolyLog[2,z] 

p n r^'^P = 

(1 — z exp(p))^ z 

= exp(p]_^dp ^ _ PolyLog[2,-z] 

(l + zexp(p))2 ^ z 



r 



The coefficients in the perturbative expansion of log(Zg) ^(3) _ ^ 21og(2) 



and log(Z/) ,are as follows; Stt^i?/ 

2^.gr^'"' ^1 ^"^'^ fc=i 
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27r2gy 


2^ 




7r3 




3 


V 


7r3 



Af^^ = — :5 for giuons and quarks 

27r^(7/ 6 



where = 8 x 2 and g/ — 6Nf are the degeneracy factors 
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